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ABSTRACT 

Context. During the last few years eclipse exoplanet spectroscopy havs yielded detection of H^O, CH4, CO2 and CO in the atmosphere 
of hot jupiters and neptunes. In the same time, ~40 likely large terrestrial planets are announced or confirmed, two of which are 
transiting, and another deemed habitable. Hence the potential for eclipse spectroscopy of terrestrial planets with the James Webb 
Space Telescope (JWST) has become an active field of study. 

Aims. We aim to extensively explore the parameter space (type of stars, planet orbital periods and types, and instruments/wavelengths) 
in terms of the signal-to-noise ratio (S/N) achievable on the detection of spectroscopic features with the JWST. We also wish to 
confront the information on the S/N to the likelihood of occurring targets. 

Methods. We use analytic formula and model data for both the astrophysical scene and the instrument, to plot S/N contour maps, 
while indicating how the S/N scales with the fixed parameters. We systematically compare stellar photon noise-only figures with 
ones including detailed instrumental and zodiacal noises. Likelihood of occurring targets is based both on model and catalog star 
population of the solar neighborhood. 

Results. The 9.6 pm ozone band is detectable (S/N = 3) with JWST, for a warm super-earth 6.7 pc away, using ~2% of the 5-year 
nominal mission time (summing observations, M4V and lighter host star for primary eclipses, M5 V for secondary). If every star 
up to this mass limit and distance were to host a habitable planet, there should be statistically ~1 eclipsing case. We also show that 
detection in transmission of the 2.05 pm CO2 feature on the 6.5 M ffi exoplanet GJ 1214 b is feasible with the Hubble Space Telescope 
(HST). For the low and the high bound of the likely atmospheric mean molecular weight, respectively just one transit or the whole 
HST yearly visibility window (107 days) is required. 

Conclusions. Investigation of systematic noises in the co-addition of 5 years worth-, tens of days separated-, hours-long observations 
is critical, complemented by dedicated characterisation of the instruments, currently in integration phase. The census of nearby 
transiting habitable planets must be complete before the beginning of science operations. 

Key words. Physical data and processes: molecular - Techniques: spectroscopic - Methods: analytical - Planets and satellites: 
atmospheres - Infrared: planetary systems - Galaxy: solar neighborhood 

1. Introduction outside the terrestrial regime. The prospect of extending the 

spectroscopy techniques above to the emerging eclipsing hab- 
Eclipsing (transit) exoplanet spectroscopy with the Hubble itable planetS; wit h the James Webb Space Telesc ope (JWST), 
Space Telescope (HST) and Spitzer has enabled the detection has been propose d dCharbonneau & Demingl2007h . The awaited 
5^ ■ of molecular signatures (H 2 0, CH 4 , CO, and CO ?i ) in the at- per f orm ance of primary and secondary ecli pse spectroscopy 



o 

(N 



O 

m 

6 



> 
OO 

(N 
O 

q 

00 
o 
o 



X 



^ w± m^uiai jigimiui^a ^'2". ^' v - v -'' ^^i! "jj performance or primary and secondary eclipse spectroscopy 

CO ■ mosphere of hot giant extrasolar planets (jTinetti et al l|200ff for habitabl e exo planets is being st udied ilBeckwithl I2008 
iGnUmairetal. 2008; Swain et al. 2009; Stevenson et al.||2010|). i Sea „ er et al 112004 ne.mW et al JbOOQl Raner et al snhmitteH - 



ISeager et al.ll2009l; iDeming et al.ll2009l Rauer et al., submitted - 

Extrasolar planets in the 1-10 Earth mass range (generally desig- rajjk) hereafter). JWST is scheduled for launch in 2014. 
nated as "super-earths") have been discov ered through Dopple r 

surveys, one of which may be habitable ( |Mayor et al. ||2009i). In M& paper we extensively explore the parame ter space 

However, past decade projects for characterisation of such plan- (steUar typeSi planet orbital period and ty and instru . 

ets (DARWIN -| Cockell et al.|2009t TPF -| Traub et al. | 2007|) are me nt/wavelength) in terms of the signal to noise ratio (S/N) 

not technologically ready yet for implementation. achievable on the detection of spectroscopic features, for pri- 

The recentl y discovered planet GJ 1214b mary or secondary eclipses, with the JWST. Because the S/N 

dCharbonneau et al.| |2009J) is the first case of close (13 pc) depends on many parameters we plot contour maps, while indi- 

transiting super-Earth, even if in this particular case the current cating ho w the S/N scales with the fixed parameters. Therefore, 

planet density estimate points to a hydrogen-rich envelope, our goa i he re is not to indicate the absolute performance of ob- 

servations, of which we do not know yet the exact conditions, but 

Send offprint requests to: A. Belu rather identify the limits and the performance gradients over the 
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parameter space. Also, some combinations of parameters will 
prove inaccessible. 

2. General target modeling 

We use model stellar parameters (mass, effective tempera- 
ture, surface gravity and luminosity) in the 0.1-1.4 M range 
dBaraffe et alj|1998l) . The stars are modeled as blackbodies. We 
have tested on the example of a 0.3 M Q (M3) star that the error 
on the S/N is below 6% for the 5-15ji/m wavelength range and 
below 15% in the 0.8-5 fim range, when using a blackbody emis- 
sion instead of a model spectrum ( lHauschildt et aljfl999l) . Only 
in the 0.6-0.8 //m range (which is not explored in this work), the 
error reaches 80%. 

The temperature of a planet is computed as a function of its 
orbital distance: 



<x is the Stefan-Boltzmann constant and F+(a) is the stel- 
lar flux at the planet's location (circular orbit of radius a is as- 
sumed). We assume a full redistribution of heat (/ = 1) for pri- 
mary eclipse observations, because we are observing the limb, 
which mixes flux from both high and low latitudes, as well from 
sunset and sunrise longitudes. For secondary eclipse, a lower re- 
distribution factor is assumed (/ = 0.75), since it is the day side 
that is observed. The Bond albedo of the planets, A, is fixed at 
0.20 for all planet types considered. 

2.1. Planet types 

We consider three planet prototypes, assumed to represent three 
large classes of planets: gas giants mainly made of H2 and He, 
icy giants (Neptune-like), and large terrestrial planets. A planet 
prototype is defined by its mass, radius, and, for primary transit 
observations, its atmospheric mean molecular weight /j. 

"Jupiters". As presented in the introduction, both primary and 
secondary eclipse spectroscopy for hot jupiters is being achieved 
today. Future instruments will give access to higher spectral 
resolutions, and cooler planets, so accordingly we consider a 
Jupiter-mass planet for the present study, with n = 2g moF 1 
(mainly H2). 

"Super-earths". Terrestrial planets can be indicatively defined 
as having an upper limit on their mass of 10 Earth masses, al- 
though in particular cases plane ts with a slight ly lower mass can 
accrete a massive gas envelope dRafikovl2006l for modeling, also 
probably the case for GJ 1214b). Therefore the prototype of a 
terrestrial planet considered in this study is a 10 Earth mass-, 2 
Earth radii-planet (hereafter "super-earth"). 

Our super-Earth prototype is considered habitable when 
found wit hin the Umits of th e circumstellar habitable zone, as 
defined bv lSelsis et al.l d2007h . For all terrestrial planets, we con- 
sider an optimistic fi = 18gmoF' (water vapor dominated at- 
mosphere of a planet in the inner portion of the habitable zone), 
instead of Earth's 28 (N2 dominated). 

1 This is the albedo the Earth would have if irradiated by a low mass 
star, because the emission maximum of the latter is shifted towards the 
infrared, where planetary molecular absorption bands are important. 
Also, a 0.1 difference in the albedo produces only a 3% difference in 
the planet's equilibrium temperature. 



"Neptunes". Neptune-mass planets represent an intermediate 
between the habitable case and the Jupiter-mass case. For this 
prototype we consider fi = 6gmol -1 (H2, He and 10% of heav- 
ier elements). 

The mean molecular weight we use is actually higher than 
that of Neptune and Uranus, which is close to 2gmoF 1 . The 
reason is that atmospheric escape is likely to deplete the amount 
of hydrogen on hot neptunes. The escape parameter for a given 
species i is X, = R p /Hi where R p is the radius of the planet 
and Hi is the individual scale height of the species i, calculated 
at the exospheric temperature. Escape becomes important for 
X < 15, while the gas is tightly bounded to the planet for 
X > 30. For Neptune and Uranus, the exospheric temperature is 
700-800 K and Xu is in the range 35-45. For warmer planets, the 
exospheric temperature of a Neptune-like planet can be much 
higher, as it roughly scales linearly with the stellar flux, until 
significant thermal io nization occurs, so up to a few thousand K 
dLammer et al.l [2003b . Exospheric temperatures above 2,000 K 
can safely be assumed for warm and hot Neptune-type planets. 
This would result in values of Xu below 15, and thus to a rapid 
escape of hydrogen. Therefore, for the hottest Neptune-type 
planets able to keep an atmosphere, the remaining atmosphere 
should be enriched in heavy elements or even consist mainly of 
heavy molecules (like N2, CO2, H2O). Therefore, our choice of 
6 g moF 1 is an average situation where only part of the hydrogen 
is left and should overestimate (respectively underestimate) ft 
for planets with a long (respectively short) orbital period. 

Table Q] summarizes the parameters of our planet prototypes. 
Note that the ~40 transiting exo-jupiters detected until now have 
a wide variety of densities, both lower and higher than that of 
Jupiter. Since we need to limit the number of parameters, we 
chose the Jupiter parameters as a middle case for the "Jupiter" 
prototype. 

2.2. Spectral signatures considered and their modeling 

S/N calculations can be based on the spectral features found in 
Solar System planets but this approach covers only a negligible 
fraction of the parameter space. They can also be based on syn- 
thetic spectra computed for specific atmospheric composition. 
In such case, the structure and composition of the atmosphere 
has to be modeled self-consistently by coupling radiative trans- 
fer, (photo)chemistry and dynamics. For the atmosphere of giant 
gaseous planets, elemental composition should not depart too 
dramatically from the stellar composition, although selective en- 
richments and depletion are expected to occur due to the separa- 
tion of condensed and gaseous phases, or to gravitational escape, 
in the protoplanetary disk and in the planet. For these planets, it 
is thus conceivable to produce grids of spectra covering a lim- 
ited number of parameters, as it is done for stars. But even in 
this case, producing such grids would imply some drastic sim- 
plifications (ID instead of 3D, equilibrium chemistry composi- 
tion instead of kinetics and photochemistry, simple cloud mod- 
els, decoupling of radiative, dynamical and chemical processes) 
and would suffer from the incompleteness of the required physi- 
cal/chemical data (spectroscopic data, kinetic rates). 

For low-mass rocky and icy planets, the situation is extraor- 
dinary more complex. Their spectral properties are determined 
by an atmosphere initially accreted as volatiles trapped in solids, 
or ices, of non-solar composition. This volatile content repre- 
sents a small fraction of the total planetary mass and is frac- 
tionated between the interior (crust, mantle), the surface oceans 
and/or ice sheets, the atmosphere and outer space through grav- 
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Table 1. Parameters of planet prototypes. 





Super-Earth 


"Neptune" 


"Jupiter" 


Radius [Earth radii] 


2 


3.85 


11 


Mass [Earth masses] 


10 


18 


317 


/i (atmospheric mean molecular mass)[g mol~'] 


18 
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itational escape (induced by impacts, exospheric heating and 
non-thermal processes). The composition of the atmosphere of 
a terrestrial planet at a given stage of its evolution is then con- 
trolled by geochemical exchanges between these different reser- 
voirs, tectonics, atmospheric escape, photochemistry, and biol- 
ogy if present. 

Therefore, the expected diversity of exoplanet atmospheres, 
and terrestrial planets in particular, covers a wide parameter 
space, and our current understanding of the origin and evolu- 
tion of planetary atmospheres provides very few constraints to 
guide us in this exploration. Although the use of detailed atmo- 
sphere models and synthetic spectra is essential, in particular to 
interpret spectral observations, it is equally important to allow 
ourselves to explore a much broader parameter space than the 
one covered today by self-consistent models. 

This is why we chose to base this study on a different, 
"model-less" approach, which is complementary to the use of 
detailed atmosphere models, which remains necessary to refine 
the actual S/N for a specific close-up in the parameter space (for 
instance RAU10, based on self-consistent habitable planet atmo- 
sphere models). 

Because of the reasons above, we chose here to examine the 
S/N of individual features of species, freeing us from any a priori 
on the atmospheric composition and structure. Moreover, what 
interests us here is not the absolute planetary signal flux, but the 
the detection of a spectral feature. 

Therefore, we model the detection by estimating the differ- 
ence of the planetary flux between two appropriately chosen 
binned channels, one measuring the continuum, and the other 
the flux in the absorption band of the feature. Of course, when a 
given (photo)spectroscopical observation comprising up to tens 
of channels will be fitted with synthetic spectra, the S/N on the 
detection of species will be much higher. With this definition, an 
S/N of 3 is a safe 3 cr detection (also see Section[3]below). 

In general, we chose to compute the S/N for a fiducial signa- 
ture defined by a given spectral resolution, and a contrast neces- 
sary for its detection. The way the signature contrast (or depth) 
is defined is described in the next section. 

However, we also wish to particularly emphasize the case 
of the habitable super-earths. As such, we consider some of the 
strongest infrared signatures of species present in the terrestrial 
atmosphere: 



CO2 feature at 4.3 fim. We measure this feature relative to a re- 
gion redward of 4 fim. Therefore, we consider in the calculation 
a mean working wavelength of 4 fim. (and an effective width of 
0.4/rni, so 7? = 10). 



CO2 feature at 1 5 fim. Since this a filter observation for JWST, 
the modeled width of the feature will be specified in the appro- 
priate section below (Section l4.1.2b . 




h=n^ 



Fig. 1. Depth of a planetary spectral feature in transmission 
(opacity height difference). is the feature's central wave- 
length, and X\ is the reference channel wavelength (continuum). 
k is the Boltzmann constant, r eq the equilibrium temperature of 
the atmosphere, fi its mean molecular weight, and g the surface 
gravity, n is the relative strength (usually n — 3, see text). 



2.3. Types of transits 

Depending of the type of transit, we use several assumptions to 
compute the planetary spectral feature depth. 



2.3.1 . Primary transit. 



We use the same formulas as iBeckwithl d2008l) for the plane- 
tary spectral feature photon count (we consider the additional 
background and instrumental noises as indicated in Eq. [2j. We 
chose the difference in atmospheric opacity height between the 
in- and out-of-band channels to be « = 3 atmospheric scale 
heights H = kT eq /fig (Figure [TJ, k being the Boltzmann con- 
stant. Consequently, the S/N scales with and n. This value 
has been observed for hot jupiters between adjacent spectral bins 
(even though larger differences in the apparent radius have been 
measured over entire spectra, see previous discussion on S/N 
definition on this page). Another way of seeing our modeling is 
as an achievable "resolution in amplitude". A "n = 3 sampling" 
should actually enable to detect opacity-radius variations of the 
planet over extended wavelength ranges (i.e. spectra) with a "bit 
depth" that could be handy if disentanglement of the signatures 
of multiple species is required). « = 3 is also a high value for 
the Earth case, where greatest opacity height difference is 4 H 
for the 15 fim CO 2 and the 9.6 //m O3 bands , and 5 H for the 
4.3 fim C0 2 band (iKaltenegger & Traub 20091) . 



O3 feature at 9.6 fim. The considered width is 0.5 fim (so R = 
20). 



2 Change of 1 order of magnitude in atmospheric pressure 
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2.3.2. Secondary transit 

Secondary eclipses allow to investigate both thermal emission 
and star reflection from the planet. We wish to have an idea of 
the weight of each phenomenon in the planetary flux, although 
they are undistinguishable in an observation. 

Emission. The radius of a transiting planet is known from the 
primary transit at short wavelengths, where the thermal emission 
from the night side of the planet can be safely neglected. The 
emitted spectrum measured at the secondary eclipse can be con- 
verted into a brightness temperature, T\,(A). The depth of spec- 
tral features in the emission spectrum depends on the vertical 
temperature profile (a isothermal atmosphere emits a featureless 
blackbody spectrum whatever its composition). At a given wave- 
length, the signal comes from the altitude za corresponding to 
an optical depth t,\ sb 1, and the associated brightness temper- 
ature is the physical temperature at this level: T\,(A) = T(za)- 
The altitude za is a function of the abundance profiles of molec- 
ular absorbers. For a given wavelength domain, the amplitude 
of the temperature variations within the altitude range spanned 
by za gives the upper limit on the depth of the observable fea- 
tures in the brightness temperature spectrum. Emission features 
can however be even stronger for a non-LTE atmosphere, where 
for instance fluorescence occur s, like it might be the case for 
HD 189733 b dSwain et alj2010h . 

The spectral signal we want to detect within the noise can 
thus be expressed as a "resolution in amplitude" (precision) to be 
achieved in measuring the brightness temperature (at the spectral 
resolution of the considered feature). We therefore consider the 
depth of the spectral feature in emission to be the difference in 
the measured planetary thermal emission between the in- and 
out- of feature channels (Figure|2]i. 

For the fiducial signature, in order not to be compelled to 
any assumption on the detailed composition and vertical struc- 
ture of the atmospheres, we chose to define the aimed resolution 
as a fraction a = 20% of the equilibrium temperature com- 
puted for the planet. Compared to the Earth case (7^ = 255 K), 
a = 20% corresponds to about the highest temperature con- 
trast (A7b = 50 K) that can be observed at low resolution on 
a disk-averaged spectrum. For this observation, 7b can reach 
near-surface temperatures of about 270 K in the 10-1 lyum at- 
mospheric window, and 220 K in the 15j im CO2 band prob- 
ing th e lower stratosphere (see for instance TChristensen & Pearl 
Il997l) . Slightly higher relative contrasts of brightness temper- 
ature have been modeled for hot and habitable Super Earths 
(RAU10) and significantly higher rela tive contrasts have actu- 
ally b een observed on hot exoplanets (M adhusudhan & Seagerl 
120101) . This value of 0.2 r eq for the brightness temperature res- 
olution appears thus as a reasonable goal to study exoplanets in 
general, and possibly the minimum required precision to search 
for atmospheric signatures habitable planets. 

We chose T eq to be the midpoint of the brightness tempera- 
ture variation, and not the upper bound (continuum temperature). 
Indeed, depending on the wavelength, atmospheric or surface 
temperatures probed by the observations can be higher or lower 
than r eq , which is a mean value. For instance, the disk averaged 
continuum emission of the Earth in the 8 and 12jum windows 
has a 7b significantly higher than T eq (due to the greenhouse ef- 
fect), except within the 9.6/iin O3 band, where 7b is lower than 
T 

A eq- 

The 4.3 A< m CO? feature is modeled in emission as having 
AT b = 100 K ( lPailletll2006l Fig. 8.9). The 15jum C0 2 and the 
O3 features are modeled with ATb respectively 60 and 30 K. 



However, around low-mass stars, the 15 fim CO2 band can be 
as deep as 100 K (RAU10). Also, O3 is the main absorbent re- 
sponsible for the stratosphere temperature inversion; some tem- 
perature profiles and O3 mixing ratios may produce a signature 
stronger or weaker than 30 K. It must be noted however (valid 
also for the primary transit case), that a habitable super-Earth 
around a low-mass star may requir e a dense CO? atmosphere in 
order to keep an atmosphere at all dJoshi et al.lll997t[Scalo et al.l 
12007b . This C0 2 may swamp the O3 signal dSelsis et alj|2002l) . 

Reflection. For reflected light spectroscopy at the secondary 
eclipse, we define our signal as a 50% difference of the specific 
(i.e. /1-dependent) planetary albedo, between the in- and out-of- 
feature channels. This is considered at 1 //m, and for a feature 
width of 0.1 jum (R = 10). 

2.3.3. Invariance of the primary transit photon count. 

Interestingly enough, we note that the fundamental physical in- 
formation for primary transit transmission spectroscopy, that is, 
the number of stellar photons traversing an atmospheric scale 
height, for a single transit, is a quantity independent from the 
planet's period P. 

Since, 7%.(a) = oT\ (R^/a) 2 , with R+ being the stellar ra- 
dius and T+ the stellar temperature, T 7 * oc P~ 4 ^, where P is the 
orbital period of the planet, therefore r eq oc P ~'/ 3 . For primary 
transits, the solid angle of the opaque annulus corresponding to 
an absorption spectroscopic feature, over a given bandwidth, is 
AD = 2R p 3H/d 2 . So AO. oc P~ i/3 . The planetary spectral sig- 
nature flux is F p (A) = AQB(A, 7 1 *), where A is the wavelength 
and B{ , ) is the Planck function. So F p (A) oc p -1 / 3 . The photon 
count for a planetary transit is n p oc F p (A) r with t the transit 
duration, r = PR+ / na, so r oc P 1 ^ . So « p is independent of P 
(in the frame of the approximations above). 

3. Signal of the spectroscopic feature and noises 

In a simple model of exoplanet eclipse observation, the plan- 
etary flux is the difference between the estimates of the in- 
transit flux and the out-of-transit one (both assumed constant 
with time). Assuming absence of correlation between the two 
(stellar photon noise dominated), the variance of the planetary 
flux is therefore the sum of the variances of the two estimates. 
Current observations of a single eclipse already have calibration 
precisions of the same order as the stellar photon noise, whether 
they cover the out-in-ou t sequence with m ultiple telescope point - 
ings dSwain et al. or a single one dGrillmair et al.ll20O7h . 

If we acquire photon counts for the transit and outside of the 
transit over the same maximum available time period (i.e the 
transit duration), these two photon counts can be considered as 
having same variance. Increasing the out-of-transit integration 
time reduces the variance of this term, in comparison with the 
in-transit one, if the only variability source is the stellar pho- 
ton noise. Stellar oscillations however contribute to the variance 
of the constant star flux estimate. This is why model fitting is 
normally used to estimate the planetary signal in transit spec- 
troscopy observations, and should be considered in a future it- 
eration of th is work. Post-detection m ethods for stellar variabil- 
ity filtering dAlapini & Aigrainll2009l) could be particularly well 
suited for massively co-added transits, around the active, lowest 
mass M dwarves (see below). 

Also, we have seen previously that we consider the differ- 
ence of the planetary flux between the in-feature- and the out-of- 
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Fig. 2. Depth of a spectral feature in emission (see text for discussion). 
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feature channels. Again, the variance on this estimate is the sum 
of the variances of the two terms (supposing uncorrelated noise 
between channels, such as photon noise). The in-feature binned 
channel width, AA, is constrained by the width of the feature (for 
instance 0.5 jim for the O3 band at 9.6 //m). The out-of-feature 
binned channel (possibly divided in two, flanking the feature), is 
chosen depending on the profile of the spectrum (optimistically 
assuming that a clean continuum can be defined). 

Following these considerations we chose a simplified model 
where the out-of-transit observation time is equal to the dura- 
tion of the transit. We also consider equal out-of-feature and 
in-feature binned channel widths. We thus compute a signal-to- 
noise-ratio on the spectral feature detection: 

g planetary spectral feature photon count for 1 transit 



V 4x (<„ 



+ erf, 



J 



where « are the photon counts for each subscript source, for the 
duration of a transit ("zodi" stands for the zodiacal light con- 
tribution, "RON" for readout noise). Exo-zodiacal light was not 
modeled. First, the contribution to the noise of an exozodiacal 
cloud (viewed in its entirety), similar to the solar one, is negli- 
gible when compared to the stellar photon noise. For systems 
with very high dust levels, its brightness variability over the 
time scale of the observation has no foreseeable source. Second, 
little is known today about the statistics of exozodiacal dust 
densities around nea rby stars, which however an ongoing effort 
(IColavita et al.ll2009l) . 

The JWST features a primary mirror of D _ 6.5 m diameter 
and a throughput before instrument of 0.88 ( iDeming et alj|2009l 
hereafter DEM09). It will be equipped with several instruments 
potentially enabling the molecular eclipse spectroscopy we are 
considering Furthermore, each instrument has different obser- 
vation modes (filter photometry, as well as low and intermediate 
resolution spectroscopy), that can be used for exoplanet spectral 
characterisation. 

3.1. NIRSpec 

In the near-infrared (N IR), we consider the NIRSpec instrument 
dRauscher et al.ll2007l) . Its performance for primary transit spec- 



3 Deming et al. Exoplanet Task Force White Paper, 2008, table at 
end. 



troscopy in the R = 1, 000 mode, for the detection of water and 
CO2 has already been studied (DEM09). We therefore focus here 
on the i? = 100 mode (0.6-5 yt/m), which could potentially yield 
a higher throughput and less readout noise, but at the expense 
of more saturation. This mode is therefore better suited for the 
faintest target stars, around which spectral characterization of 
super-earths will be most efficient. We actually show below that 
the saturation will not be a limiting factor, considering the num- 
ber of target stars of a given type within a given distance from 
the Sun. 

The overall throughput (including quantum efficiency) and 
the resolution function were provided by P. Ferruit (CRAL, 
Lyon, France). Since exoplanet transit observations will be done 
in a pseudo slitless mode (custom 1.6" 2 opening), we scale the 
provided resolution function (which was computed for the 0.2" 
slit) with the point-spread-function (PSF) /l-dependent size. For 
wavelengths under 1 fim we use the PSF size at this upper bound, 
to account for distortions that become non negligible below this 
value. The resulting resolution curve serves to compute the num- 
ber of pixels in the spectral feature channel, hence the readout 
noise. The pixel scale is 0.1", the readout noise is 10e~/pixrms, 
and the well capacity is 60,000 e . Since the PSF is undersam- 
pled at the shortest wavelengths, we suppose that the spectrum 
axis is centered between two rows of pixels. A defocus mech- 
anism could however be present to mitigate the undersampling 
problem. 

In order to compute the brightest pixel, we use the resolution 
function at the shortest wavelength (maximum star emission), by 
assuming that there is no diffraction in the dispersion direction. 
This is an accurate model in terms of energy distribution over 
short dispersions which are part of a larger spectrum (the reso- 
lution function is computed for a 2.2 pixel size of the resolution 
element). In the spatial direction, we use a simple triangle model 
of the center of our dispersion-anamorphosed PSF (83% of the 
total energy). The maximum number of electrons in the brightest 
pixel is used to compute the readout rate. 

To determine the read time, we compute the length of the 
NIRSpec R = 100 spectrum from the resolution function 
(356.2 pixels), rounded to the upper power of 2 (512). The width 
(in the spatial direction) on the detector is 2 xA mRX /D (A max = 
5 /im here), rounded in pixels to the closest upper power of 2 
(8 pixels). We assum e the read mode is MULTIACCUM-2xl 
(Rauscher et a 1120071) . meaning that we only have to account for 
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Fig. 3. Brightest pixel saturation time for the full R = 100 
NIRSpec mode, as function of star mass (star at 10 pc). Dotted 
level: spectrum window reset time. 



the reset frame time (the detector is read non-destructively up- 
the-ramp). We reduce the effective photon collection time over 
the transit by this amount (the remaining fraction is called duty 
cycle). 

Figure [3] shows the brightest pixel saturation time as a func- 
tion of stellar type, and the full spectrum window reset time 
level, assuming each pixel requires 10/is to reset. Smaller de- 
tector windows can be defined (possibly alternately observed), 
depending on the sought spectral signatures. However, we have 
not modeled such observations, so we limit our NIRSpec plots 
at the bound defined by the intersection of the reset time and the 
saturation time: 1.4 M Q at lOpc (four F4 stars only: Sirius A, 
Altair, Formalhaut, Vega and Procyon A). The limit is 1.05 M Q 
at 5 pc (three stars only: a Cen, Sirius A and Procyon A). 

Finally, the PSF being undersampled at the shortest wave- 
lengths, the spectrum's jitter over the pixel grid (7 mas rms, 
DEM09) is likely to constitute the major noise source in the 
NIR; however, this modeling is beyond the scope of this work. 
For instance, mapping all detector pixels in the dispersed custom 
window, before flight, should enable processing techniques that 
reduce this noise contribution. 

3.2. MIRI 

The mid-infrared ( MIR) range is covere d by the Mid-InfraRed 
Instrument (MIRI, IWripht et all l2004l) . Its filter photometry 
mode performance has been already studied for the detection 
of the CO2 band at 15/im, for secondary transit spectroscopy 
(DEM09). However, the current filter set is not optimal for the 
detection of the 9.6yum O3 band (the JM_3 imaging filter cov- 
ers 9-11 //m, four times the width of the feature). We therefore 
consider here the 5-11 fim, low resolution (R = 100) spectrome- 
ter mode (LRS). The optical transmittance is assumed to be 0.4 
for both modes. Both modes use the same detector with a pixel 
scale of 0.11", a readout noise of 19e~/pix rms, and a well ca- 
pacity of 10 5 e~. Again, we do not consider the effect of the 
aforementioned instrumental jitter on MIRI (PSF undersampled 
shortward of 7 /an), although it was checked in filter mode on 
one simulation that its effect is negligible compared to other lim- 
itations (Cavarroc et al., forthcoming). 

We assume a total spectrum length of 194 pixels in the dis- 
persion direction (wavelength pixel registering provided by S. 



Ronayette, CEA). We make the assumption that the read mode 
is similar to that of NIRSpec (above). 

The extension of the feature on the spectrum is computed 
by a simple proportionality between the width of the spectral 
feature and the total length of the spectrum. No saturation 
occurs for MIRI at 10 pc over the considered 0.1-1.4M o range 
of stars. 

For both instruments, a 0.03 e~ s _1 pix~' dark current noise is 
considered. Uniform background noise sources are calculated 
using the pixel scale of the detectors. For the instrument's 
thermal emission we use a temperature of 45 K and a 0.15 
global emissivity. 

3.3. Zodiacal light 

We use an implementation by R. den Hartog of a paramet- 
ric model by O. L ay of the Kelsall local zodiacal cloud model 
dden Hartogl2005l) . The most pessimistic ecliptic latitude (J5 = 0) 
is used, but in the anti-solar direction, which represents a good 
average of the all-sky distribution. 

4. Results 

We now present the computed S/N for the different types of tran- 
sits, for different types of planets, and for different wavelength 
ranges. For the smallest planets, achieving a significant S/N will 
require cumulating data from multiple transits, in which case the 
S/N scales with the square root of the number of transits, pro- 
vided that (instrumental) noises are not correlated between the 
successive transit observations. Hence, the detection limit is ul- 
timately the mission life-time (5 years), which limits the S/N 
that can be achieved on longer period planets. It must be empha- 
sized that such an hypothetical observation, while having a total 
observation tim e only a magnitude o ver the longest exposures 
made until now dBeckwith et al. 200©, suffers from the risks in- 
herent of being distributed over a 100 times larger duration. In 
other words, if JWST becomes inoperable after 2.5 years of op- 
erations, the 1 % of the mission time dedicated to acquiring data 
on the planet would be lost, since, the data would yield an insuf- 
ficient S/N. 

It must be noted that the yearly target visibility for JWST 
reaches 100% only for targets with ecliptic latitude higher than 
85°. The visibility is lower than 100 days per year for latitudes 
up to 4513- We find the yearly mean sky visibility to be 149 days. 
For all following super-earths plots, we therefore multiply the 
number of transits occurring over the 5-year mission time by 
the corresponding fraction, and obtain the effective number of 
observable transits. For all planet types, we also ponder for the 
unknown impact parameter of the transit, by further multiplying 
the transit duration (equatorial) by n/4 (~0.79). 

Unless mentioned otherwise, all our examples here are cal- 
culated for a system at a distance d = 10 pc, the signal to noise 
scaling linearly with the inverse of the distance. Jupiter-mass 
planets are studied at 50pc (see Section|6]i. 

4. 1 . Primary transit 

The main interest of primary transit spectroscopy is that it yields 
a spectrum even in the case of an atmosphere which is isother- 
mal, or which has a low temperature gradient (this would not be 

4 http://www.stsci.edu/jwst/overview/design/field_of_regard.html 
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case if this atmosphere where observed in emission, during a sec- 
ondary transit). Also, the primary transit signal is proportional to 
T eq , whereas this is true in emission only in the Rayleigh-Jeans 
regime. As a consequence, some (A, r eq ) couples (cold planets 
at short wavelengths) may have an undetectable emission flux, 
but could be characterized in primary transit (see below). 

4.1.1. NIRSpec 

Figure|4]shows the S/N for the detection of the 4.3 /im CO2 band 
as modeled in Section |2~2"1 for a super-Earth (this figure serves 
as a template for following similar plots; refer to its caption for 
conventions). All primary transits available on average (see be- 
ginning of current section) are cumulated over the 5 year mission 
time, with stellar noise only (left) and the instrumental and zo- 
diacal noises (right). We have checked on several examples that 
our figures with stellar noise only compare well with RAU10. 
The main noise influence here is the readout-noise. CO2 can 
be detected on super-earths in transmission. The distance to the 
star being fixed in this type of plot, the decrease of S/N towards 
brighter (heavier) stars when considering instrumental noises is 
due to the increasing weight of the reset time in the duration of 
each exposure, hence the reduction of duty cycle. We limit our 
study to the main sequence dwarves M through F. While the lat- 
ter, more extended stars should provide an increased duration 
for collecting photons (transit duration), the duty cycle displays 
a maximum at 0.6 M . The position of this maximum is 0.5 M Q 
(MO) at 5 pc. 

In order to explore the sensitivity of the instrument at var- 
ious wavelengths (for other potential spectral signatures), we 
consider the fiducial spectral signature (« = 3 scale heights H, 
and of constant spectral resolution of 20 over the wavelength 
range - Figure[5]l. The S/N scales linearly with the inverse of the 
square root of the resolution. The dominating noise source is the 
readout noise, set by the saturation time of the brightest pixel of 
the spectrum. The variation with stellar type of the wavelength 
of maximum emission is clearly visible. The abrupt reduction 
of NIRSpec throughput towards the shorter wavelengths is also 
visible when comparing the star-only-noise plot (left) with the 
full modeled noise one (right). 

Figure [6] shows that Neptune-size planets will require sev- 
eral transits to achieve spectroscopy at the maximum resolu- 
tion of the mode (R = 100). For jupiters, the performance on 
their S/N at 50 pc (Figure [7]) is similar to that of the neptunes at 
lOpc above. However, given the greater distance, the saturation- 
induced curbing of the S/N with star mass is shifted towards 
brighter stars: therefore Jupiter-size planets can be characterized 
(at R = 100) with only one transit up to a period of a few hun- 
dred days. 

4.1.2. MIRI 

Figure [8] shows the S/N for the detection of the 9.6 ji/m O3 fea- 
ture in our super-Earth prototype, for observation of all the pri- 
mary transits available on average over the 5 year mission time, 
with stellar noise only (left) and all the above MIRI and zodiacal 
noises (right). As seen above for the NIR spectroscopy, habit- 
able planets can be characterized only around low mass stars. 
Additionally, because the O3 feature is difficult to detect, we 
calculate this and all subsequent plots for this species at 6.7 pc 
(value derived from occurring statistics for transiting habitable 
planets, see Section|6]l. It can be seen that the ozone feature will 



be detectable in transmission only for warm habitable planets 
around the lowest mass M dwarves. 

It should be noted that the chosen value for (lSgmoF 1 ) 
implies an atmosphere dominated by water- vapor, corresponding 
to a super-Earth with a sufficient water reservoir near the inner 
edge of the habitable zone. O3 detection as well as its interpreta- 
tion in terms of biosignature are problematic within a JrbO-rich 
atmosphere: 

- photochemical productions of H, OH, HO v associated with 
H2O photolysis drastic ally limit the build-up of an ozone 
layer dSelsis et al.ll2002l) , 

- collision broadening far-wi ng absorption by H 2O can screen 
the O3 signature (see Fig. 4. ISelsis et aT1l2008l) . 

- the presence of H2O above the tropopause yields enhanced 
loss of hydrogen to space and the abiotic build-up of oxygen 
left overs, making the indirect biological origin of O3 doubt- 
ful (ISelsis et al.ll2007h . 



In order to compare with forthcoming work from Cavarroc 
et al., we also consider the CO2 feature at 15yum. We use the 
filter iM_6, 3yum wide, centered at 15yum. No saturation occurs 
for this observation for the range of considered stellar masses, 
assuming that we can read only the PSF window (2 times the 
size of the PSF corresponding to the upper wavelength bound of 
the redmost filter (7JVL6). For the out-of-feature reference, we 
use the IMA filter, running from 10.95 to 11.65yum, that is 3 
times narrower than IM_6. We therefore substitute the factor 4 
in Eq.|2]by ( 3 /g + l) + 2 ~ 3.33, but only for the bandwidth- 
dependent terms (i.e stellar, zodi and thermal). The effective 
bandwidth considered is therefore that of the out-of-feature filter 
{IMA, 0.7/mi), and the effective wavelength is 11+15 /2 = 13 /mi. 
Figure [9] shows the S/N for this feature. Given the uncertainties 
in our model, we can see that CO2 at 15 /mi will be difficult to 
detect in transmission even for warm habitable planets around 
the lowest mass stars. 

At constant stellar mass, the reduction in S/N is a compos- 
ite of a) the uniform effect of the readout noise, b) the reduced 
quantum efficiency of the MIRI detector for the shorter wave- 
lengths, c) the thermal and the local zodiacal contribution to- 
wards longer wavelengths (FigurefTOli. It can be seen that perfor- 
mance increases towards shorter wavelengths, where additional 
signatures such as methane and water bands are situated. 

Thanks to their lower density^, because of their high hy- 
drogen content, Neptune-size planets can be spectroscopically 
characterized from a single or a couple of transit observations. 
Figure [TT] shows the S/N for a single transit of our prototype 
Neptune planet, towards the most unfavorable end of the MIRFs 
LSR wavelength range. The S/N for our Jupiter prototype in the 
(wavelength-planet period) space, around a solar type star, and 
at the full resolution of the instrument, is presented in Figure[T2l 

It can be noted that for the giant planets, for NIRSpec as 
well as for MIRI, the dynamic of the S/N over the parameter 
space is not very strong, implying that these planets can be 
characterized in a wide variety of cases with a fairly constant 
number of cumulated transits. 



5 Thus both high radius and high H. 
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Fig. 4. S/N on the detection of the 4.3yum CO2 signature in a super-Earth at lOpc, for observations of all the primary transits 
available on average over the 5 year mission time. For this and all following S/N plots: left plot is with stellar noise only and right 
plot is with instrumental (here NIRSpec) and zodiacal noises. For this plot and following cases of observations over the whole 
mission time (super-earths): a) the habitable zone is plotted in gray, b) the fraction of the mission time, accounting for 2 transit 
durations ("1/2 out + in + 1/2 out-of-transit"), plus the fixed 65 min JWST slew time, is on the right axis. 
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Fig. 5. Same as Figure |4] but in the (wavelength - star mass) parameter space. The planet's distance to the star is such that it 
receives the same amount of energy as Earth (1 AU from the Sun). The width of the fiducial spectral feature used in the computation 
corresponds at each wavelength to a fixed resolution (R = 20). 



4.2. Secondary transit - emission 

It is interesting to note that the atmospheric species detectable 
in an emission spectrum are, by definition, greenhouse gases, 
which affect the planet's climate. 

4.2.1. NIRSpec 

Even when integrating over the whole mission time, the 4.3 fim 
CO2 feature is not detectable in emission on habitable super- 
earths, even those on the inner edge of the habitable zone 
(Figure IT3b . However, the CO2 signature can be detected in 
emission on hot super-Earths around low mass stars. Let us sup- 
pose for instance a high gravity planet and a dense atmosphere, 
that would lower the scale height and prevent detection in trans- 
mission. D etection of CO2 through emission would be an indi- 
cator of the presence of an atmosphere, which is a question of 
debate for these low mass objects at short distances from active 
stars. 



Unlike the case of transmission spectroscopy, Neptune-mass 
planets cannot be characterized in emission with NIRSpec with 
only one or a couple of transits (plot not shown). It is required to 
go up to the Jupiter-mass scale to recover this ability, and only 
for the hot planets. The performance on the S/N drops drasti- 
cally shortward 2 pm (Figure [T4l. For these cases, the strength 
in emission of our R = 100-equivalent wide fiducial spectral fea- 
ture is A7b = a T eq with an optimistic a = 0.2. 



4.2.2. MIRI 

Figure [15] shows the secondary transit-emission spectroscopy, 
for our habitable planet prototype, for the O3 feature. With our 
current set of parameters, the signature appears less detectable 
than in primary transit (Figure [8j. 

We examine then the 15 /jm CO2 signature in emission, with 
the same modeling as in Section 14.1.21 Figure [16] shows that 
(given the considered parameters) the detection of CO2 in emis- 
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Fig. 6. S/N for a Neptune-size planet at lOpc, at the mid wavelength (3/im) and maximum average resolution of the considered 
NIRSpec mode (/?= 100), and for a single primary transit. For this figure and all following gas giant planet plots: note the 
extended orbital period scale (where applicable). 




Fig. 7. Same as Figure[6]for a Jupiter-size planet at 50 pc. 



sion at 15 //m is a little less efficient than in transmission at 
4.3 /mi (FigureHJ. 

As before, for Neptune and Jupiter-mass planets we consider 
a fiducial spectral feature, with the same a = 0.2 depth. Results 
are shown in Figure [T7] For jupiters for instance, the S/N is 
not so sensitive with wavelength as in the primary transit case 
(Figure fT2l. 

4.3. Secondary transit - reflection 

Reflection spectroscopy is considered only with NIRSpec, be- 
cause reflection is swamped by thermal emission in the mid- 
infrared. Reflection spectroscopy of super-earths is not achiev- 
able even over the full mission time, and is not shown here. Also, 
reflection spectroscopy of Jupiter-mass or Neptune-mass planets 
is achievable in one or a couple of transits only for the most il- 
luminated planets (Figure[l8]for the Jupiter-mass case). 

5. Test on a HST observation 

In order to test our model, we implemented the corresponding 
characteristics of the Hubble Space Telescope (HST) NICMOS 
instrument for observations wi th the G206 grism (slitless). We 
simulate the lSwain et "all d2008l) observation of HD 189733 b. We 



choose a feature of 2 scale heights (1,400 K assumed for plane- 
tary temperature) and one channel wide (R = 40) at 2.05 /mi. In 
this paper, we have assumed that the individual pixel responses 
are well characterized, enabling to observe very near to the sat- 
uration limit. With this setting, our HST duty cycle is very near 
100%, so we adjust the S/N by the reported 18% duty cycle of 
this observation. Our model S/N is 3.6 stronger than the one of 
the observations, which can account for the pointing oscillations 
and other optical state variations that we do not model. 

For prediction purposes, we have investigated the S/N for the 
observation of the primary transit spectroscopy GJ 1214b at the 
same wavelength and resolution (illustrative of the detection of 
the CO2 band, 3 scale heights strong). The largest HST yearly 
visibility window for GJ 1214 is 107 days. We consider a planet 
period of 1.6 days, a transit duration of 48 min, and a HST orbit 
period of 96.5 min. We assume that the target is visible only 40% 
of each orbit. We compute there are on average 59.6 (cr=2.6) out 
of the 68 within-the-window transits (88%), which are at least 
partially covered by the visible period of an HST orbiQ A to- 



6 GJ1214 has no Continuous Viewing Zone (CVZ) for HST. 
Because of atmospheric drag of the low orbit- HST, it is impossible 
to know even a couple of weeks beforehand the exact position of the 
telescope. 




Fig. 9. Same as Figure[8] but for the 15 /mi CO2 feature, and at 10 pc. 



tal of 21.3 (<x=0.5) hours of transit are covered by effective orbit 
portions, and can therefore be cumulated. The required telescope 
time is at least the double because of the out-of-transit observa- 
tions for the stellar baseline determination. Applying the above 
3.6 scaling factor to the result of our model, we obtain an S/N of 
1.8for/f = 18 g moF 1 , and > 5 for// = 6. If we consider// = 2 
(consistent with the hydrogen-dominated atmosphere currently 
assumed to explain the observed radius of the planet), S/N = 
2.5 with only one transit. The observation of only one secondary 
transit with JWST-MIRI (a = 0.2, R = 20 at 10pm) would 
yield a S/N of 3.5. 

6. Likelihood of occurring targets 

6. 1 . Habitable planets 

We have seen that transiting habitable super-earths could be 
characterized at low-resolution with a significant S/N only in 
the most nearby systems (<10pc) and when hosted by a M 
star (<O.2M ) (Figures S [3 [J3j [B] and [16]). These two cri- 
teria match well since the stellar population is dominated by 
M stars; however, we are dealing with a small number of stars 
(~300 within 10 pc). We thus calculate the occurrence likelihood 
of transiting habitable planets around M stars (Figure [T9b. We 



make use of a complete census within 6.7 pc made by RECONS 
(Research Consortium on Nearby Stars fl- For each star, we com- 
pute the transit probability for a planet receiving the same inso- 
lation as the Earth, and we can thus calculate the statistic number 
of transiting planets as a function of the distance to the Sun. This 
value has to be multiplied by 77®, the mean number of habitable 
planets per star. Within 10 pc and for M0-M9 stars, 4.6 x t] b tran- 
siting habitable planets are expected. We have excluded K, G and 
F stars because of the low S/N for habitable planets. Including K 
stars would increase the statistical number of transits by 10%. 

One should note also that rj m here is not the fraction of stars 
having a habitable planet because one star may host several hab- 
itable planets, which was maybe the case for the Sun 4 Gyrs 
ago, when Venus, the Earth and Mars were potentially habitable 
dSelsis et al.ll2007l) . The case of 775, > 1 is thus not to be dis- 
carded. 

Because the Earth is located near the inner edge of the habit- 
able zone, using its insolation to compute the transit probability 
yields optimistic numbers: if we assume a uniform distribution 
of planets with orbital distance, most habitable planets have a 
transit probability lower than the one we use. Climate models 
predict that the outer boundary for the Sun's habitable zone is 

7 |http://www.chara. gsu.ed u/RECONS 
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Fig. 10. Primary transit full spectroscopy S/N for a super-Earth at lOpc with MIRI, in the (wavelength - star mass) parameter space. 
The planet's distance to the star is such that it receives the same amount of energy as Earth (1 AU from the Sun). The width of the 
fiducial spectral feature used in the computation corresponds at each wavelength to a fixed resolution {R=20). 
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Fig. 11. S/N for a Neptune-size planet, for a single primary transit, at IQfim (MIRI), and at the full resolution (R = 100) of the 
instrument. 



located around 2 AU. Using 2 AU instead of 1 AU as a reference 
yields transit probabilities decreased by a factor 2. Therefore, 
for j] @ = 0.3 the statistical number of transiting habitable plan- 
ets reaches 1 somewhere between 8.5 and 12pc for the whole 
M0-M9 range and between 10 and 16 pc for the M5-M9 range. 



We should also note that our S/N calculations for habitable 
planets assume a 2R B radius, which may be significantly larger 
than the average value. For primary transits, the atmospheric sig- 
nal varies as /S ,s for rocky planets (because density varies with 
the mass). The S/N scales as R^ for secondary eclipses. For this 
reason also, the values we chose seem to represent an optimistic 
situation where large planets dominate the population of habit- 
able planets. 



Last, the computed S/N for a habitable super-Earths assumes 
an observational program regularly spread throughout the JWST 
mission time, and therefore a transiting target known at the be- 
ginning of the mission. 



6.2. Jupiters. Neptunes and hot super-earths 

Hot jupiters (> 0.1 jupiter masses and period < 10 days) are 
found by radial velocity ( RV) surveys around ab out 2% of F- 
G-K stars (see for instance [Cumming et al. 1 120081) . Their period 
distribution peaks at ~3 days. With these two values we com- 
pute the transit probability for individual nearby stars of a syn- 
thetic population of nearb y stars generated with the Besancon 
model (iRobin et al.ll2003ln We find that the statistical number 
of transiting hot jupiters reaches 1 at 23 pc. The closest known 
hot jupiter is HD 189733b, and is found at 20 pc. We should ex- 
pect about 10 transiting hot jupiter within 50 pc, among which 2 
are already known (HD 209458 b and HAT-P-1 1 b). This is why 
we choose this distance for our S/N contours for hot jupiters. 

Lower limits on the freque ncy of longer-period g iant plan- 
ets can be estimated from RV dCumming et alj|2008l) . A clear 
trend shows that long period jupiters are more frequent than hot 
ones, with at least 10% of the stars hosting gas giants with pe- 
riod smaller than Jupiter's one. This increase of frequency does 
not however compensate for the linear decrease of the geomet- 



http ://model . obs -besancon. f r 
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Fig. 12. Same as Figure[TT] but for a Jupiter-size planet at 50 pc around a solar-mass star, in the wavelength-planet period space. 
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Fig. 13. S/N in emission of the 4.3 /an CO2 signature in a super-Earth at 10 pc, for observations of all secondary transits available 
on average over the 5 year mission time, with stellar noise only (left) and modeled NIRSpec and zodiacal noises (right). 



ric transit probability with the orbital distance. This means that 
the number of transiting cases observed within a certain distance 
decreases as the orbital period increases. 

This conclusion may not be true for M stars, for which no hot 
Jupiter have been found but long period gas giant are detected by 
RV and microlensing. Microlensing results suggest that 40+20% 
of the lenses host a mass ive planet at the pro bed orbital distances 
(beyond the snow line, iGould et al.l [20 1 0t> . The population of 
lenses is dominated by stars in the 0.3-0.7 M range. Let 77 be 
the mean number of planets of considered type per considered 
star mass range. If we extrapolate the above result and assume 
a value 40% for 77 at the snow line0 for stars down to 0. 1 M Q , 
then the statistical number of transits reaches 1 at only 1 1 pc. 
However, despite this potentially high number of nearby transit- 
ing long-period planets, these systems are hard to find with the 
current methods, and may not be unveiled in time for JWST. 

Short period planets (<50 days) in the 5-20 Earth mass range 
have been fou nd by the HARPS RV survey around 30+ 10 of the 
G and K stars dMavor et al.ll2009bl) . Taking into account the un- 
certainty on r] and the fact that the dependency of rj upon the 
period is yet to be determined, the statistical number of transit- 
ing systems reaches 1 between 8 and 18pc for G-K stars only, 



Taken at 2.7 AU for the Sun and scaled with luminosity 



and possibly as close as 5 pc if extrapolated to M dwarfs. These 
hot low-mass planets represent promising targets for JWST, as 
determining if hot super-earths can sustain an atmosphere, and 
of what composition, is a key scientific question. 

7. S/N scaling with various parameters 

Table [2] summarizes how the S/N scales with the different pa- 
rameters considered constant in the contour maps. Each param- 
eter may scale different components of the S/N (for instance, the 
distance to the star d scales the stellar photon noise S/N). 

8. Conclusion 

In this paper we have computed S/N for the detection of spec- 
tral features in exoplanetary atmospheres through eclipse spec- 
troscopy with JWST. We insist it is a S/N on the detection of 
spectral feature, and not on the measured value of the flux. The 
S/N is represented as function of exoplanet parameters (size, in- 
solation, host star, transit duration and frequency) and observa- 
tional ones (resolution, wavelength). Our spectral features are 
modeled by a couple of parameters only (strength, width), in or- 
der to better explore the parameter space, and identify the re- 
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Fig. 14. Same as Figure [13] but for a single secondary transit of a Jupiter-size planet, at 50 pc, around a solar type star, in the 
wavelength-planet period space, and at the full resolution of the NIRSpec R — 100 mode. The depth of the fiducial spectral 
signature is a T eq with a = 0.2. 
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Fig. 15. S/N in emission for the 9.6 /urn O3 signature, with the MIRI instrument for a star situated at 6.7 pc. 
Table 2. S/N scaling with various parameters 



Parameter 




Scales as 


Distance to star 


d 




Collector diameter 


D 


D 


Number of observed transits (uncorrelated noise) 


N 


Viv 


Planet radius (super-Earths) 


R P 


.ftp 8 (primary) or (secondary) 


Atmospheric mean molecular mass 


I* 


7„ 


Number of scale heights (feature opacity) 


11 




Resolution 


R 





gions of interest where detailed atmosphere models and syn- 
thetic spectra can bring further insight. We systematically com- 
pare a stellar photon noise-only scenario with one containing 
background and instrumental noises. 

One primary transit observation with JWST-NIRSpec 
(R-100 mode) will enable to attain S/N = 3 at 3jum on giant 
planets. At 50pc, this result can be accomplished for jupiters 
with periods up to ~300 days around G and F stars (would they 
be detected), and up to 30 days around K stars. For neptunes 
(lOpc away) the same S/N (same wavelength and resolution) is 



obtained for planets with periods up to again ~ 300 days but 
around K and G stars. 

In the MIR (MIRI), at 10/im, S/N= 3 (R = 100) can be 
obtained by summing <10 transits, for neptunes up to 100 days 
period. For jupiter-mass planets the period limit for MIR S/N = 4 
with 4 transits decreases from 30 days at 5 jum to 4 days at 1 1 //m 
(solar host assumed). In secondary eclipse MIR observations, the 
limit on the contrary increases with wavelength, from 20 days 
at 5yum to 80 days at 11 /im for jupiter planets. For neptunes, 
the limit is roughly uniform with wavelength: ~ 15 days. In the 
NIR, neptunes become difficult to characterize even by summing 
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Fig. 17. S/N for a single secondary transit of a Neptune-size planet at lOpc (top row) and a Jupiter-size one at 50 pc (bottom row), 
around a solar-mass star, in the wavelength-planet period space, and at the full resolution of the MIRI low resolution spectrometer 
(R = 100). 



transits, and the jupiter planet limit (S/N = 5 with 1 transit, R = 
100) is around 4 days and only for wavelengths over 2 fim. 

In summary, JWST will better characterize exoplanets 
already characterized today in eclipse spectroscopy, but, more 
importantly, will enable to characterize at the level available 



today (or better) objects more distant, thus more numerous, 
therefore contributing strongly to comparative planetology. 

We devote particular attention to the prospect of character- 
izing habitable planets. Since achieving a significant S/N 
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Fig. 19. Number of transiting systems receiving Earth's inso- 
lation, as a function of distance and mean number of habit- 
able planets per star (t] @ ). The doted curves are based on the 
RECONS (Research Consortium on Nearby Stars) list. The solid 
curves results from rejecting stars separated from a companion 
by less then 5", as well as well as GJ581 (already surveyed for 
transits) and GL 876 (eccentric planet in the habitable zone). The 
gray profiles are a d 3 fit (±cr) to the list with close binaries re- 
moved. 



forts on the most nearby stars, whatever the brightness. We also 
wish to stress that the census transiting habitable planets around 
M dwarves should be as complete as possible by the beginning 
of JWST operations. Otherwise, we may be faced with the sce- 
nario of choosing to start an observational program worth 2%0 
of the 5 year-mission time, only to find after 2.5 years (halfway) 
that there is a more interesting target. 

While a timely precursor to dedicated observatories as 
Darwin or TPF-I/C, habitable exoplanet eclipse spectroscopy 
with JWST can not reach the science objectives of the former. 
K and G stars are out of reach , and M star planets have a lot of 
habitability issues IScalo et all d2007l) . Photon collection rate is 
limited by duration and frequency of the transit, and (not undis- 
sociated from the previous) the technique can not give access to 
enough targets to provide comparative planetology statistics. 

Further work will include considering the exact variances of 
the light curve fitting techniques actually used for the exoplanet 
spectrum calculation (ICarter et al.ll2008l ). We have not modeled 
limb darkening (for primary transit), star spots, and stellar vari- 
ability which is likely to be peculiar for M dwarves. We rec- 
ommend detailed preflight characterisation of the detectors. For 
one, this would enable to mitigate for pointing oscillations. Also, 
we have assumed here that we work near the saturation limit of 
the detectors, where their behavior is highly non-linear. Our code 
is available on request (contact first author). 



requires to sum a large number of transits we systematically 
compute the cost (fraction of the JWST mission time) of such 
observations. Detection of O3 at S/N > 3 is feasible around 
M4-M9 stars in primary and around M5-M9 in secondary 
eclipse for warm super-earths 6.7 pc away or closer, with ~2% 
of the 5 year JWST mission time. We compute that if every 
M star out to 6.7 pc where to have one habitable planet, there 
should be ~1 transiting case. CO2 (15 et 4.3/im features) are 
also detectable (only the 4.3 yum is not detectable with secondary 
transits). Shorter (longer) wavelengths are naturally more suited 
for primary (secondary) eclipse observations. 

We showed that low mass stars that yield a high S/N on the 
spectroscopic observation of their habitable planets are also the 
peak in the number of likely transiting cases. Because of the need 
for radial velocity confirmation, most transit searches survey the 
brightest stars (magnitude limited samples), while if the purpose 
is the characterization of atmospheres, we should put more ef- 
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